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ABSTRACT 


This  report  is  the  first  quarterly  technical  report  under  this  study  of  reentry  effects 
associated  with  the  transmission  of  electronic  countermeasure  (ECM)  signals.  Illustrative 
results  for  the  aerothermochemical  properties  are  presented  for  typical  slender-body  re¬ 
entry  conditions.  These  results  were  obtained  from  an  approximate  nonequilibrium  inviscid 
shock  layer  flow  field  solution,  which  utilizes  a  combination  of  reacting  streamtube  and 
equilibrium  flow  field  solutions.  A  discussion  of  the  calculation  method  is  included. 

Literature  surveys  in  the  area  of  theory  of  antennas  and  propagation  are  summarized. 

(Types  of  ECM  antennas  in  common  use  are  discussed  in  Volume  n. )  Methods  of  calculating 
near -field  distributions  of  antennas  are  discussed  in  a  general  way.  Theoretical  analyses  of 
linear  plasma  effects  are  surveyed.  Antenna  breakdown  is  considered  in  terms  of  both  the 
static  air  environment  and  the  reentry  environment.  Hie  theories  of  nonlinear  propagation 
effects  are  summarized  and  finally  a  review  of  some  alleviation  techniques  is  given. 

In  general,  it  appears  that  nearly  all  of  the  important  reentry  effects  on  antennas  and 
propagation  can  be  reasonably  well  characterized  using  existing  theory.  However,  parametric 
evaluation  using  these  theories  will  have  to  be  completed  before  this  conclusion  can  be 
affirmed. 
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SECTION  1 
INTRODUCTION 


irv.'-' 


The  study  under  which  this  report  is  written  covers  two  main  theoretical  disciplines:  fluid 
mechanics  and  electromagnetic  wave  propagation  in  plasmas.  Generally  speaking,  when 
theoretical  calculations  are  to  be  performed,  it  is  necessary  to  obtain  the  fluid  mechanic 
effects^  upon  which  the  reentry  plasma  properties  depend  before  proceeding  to  the  electromagnetic 
propagation  analysis.  The  scope  of  work  in  both  areas  covers  literature  reviews  as  well  as 
parametric  calculations.  But  in  order  to  allow  lead  time  to  the  propagation  calculations,  the 
parametric  flow  field  evaluation  was  begun  immediately,  llie  final  report  will  cover  the 
fluid  mechanic  state  of  the  art  evaluation.  This  report  covers  the  results  of  inviscid  flow  field 
calculations  by  the  streamtube  method,  wtiicMs  described  in  detail  below. 


The  problem  of  the  propagation  of  signals  from  a  radio  wave  source  aboard  a  vehicle  re¬ 
entering  the  atmosphere  at  hypersonic  speeds  involves  many  areas  of  electromagnetic  theory 

* 

and  plasma  physics.  Antennas  must  be  considered  in  terms  of  electrical  characteristics  and 
near  field  distributions.  The  linear  effects  of  interaction  of  the  antenna  with  the  reentry 
plasma  must  be  considered  with  respect  to  antenna  impedance  and  pattern  distortion. 

Antenna  breakdown  is  an  important  consideration,  especially  at  high  transmitted  power  levels. 
This  problem  must  first  be  considered  in  the  static  environment  of  undisturbed  cool  air. 

Then  the  effects  of  reentry,  where  experimental  data  are  much  less  plentiful,  must  be 

I 

considered.  Hopefully  an  extension  of  the  well  substantiated  theory  for  the  static  environment 
can  be  carried  out  in  such  a  way  as  to  cover  reentry  effects. 

I 

Nonlinear  propagation  effects  must  be  considered  in  terms  of  attenuation  and  distortion  as 

| 

1 

functions  of  input  power  level.  Finally,  the  possible  techniques  of  alleviating  all  of  the  above  ' 
effects  must  be  evaluated.  All  of  these  areas  of  antennas  and  propagation  are  summarized  in 
the  following  sections  in  terms  of  the  significant  literature. 
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SECTION  2 

FLOW  FIELD  ANAI  YSIS 


Typical  results  from  the  inviscid  flow  fields  which  have  been  generated  as  inputs  to  the 
parametric  study  of  the  electromagnetic  wave  propagation  problems  are  presented  and 
briefly  discussed.  In  both  the  selection  of  these  results  and  the  accompanying  discussions, 
emphasis  has  been  placed  on  the  nose  bluntness  and  altitude  effects  on  the  nonequilibrium 
flow  properties,  and  particular  attention  has  been  directed  toward  the  contrast  between 
chemical  equilibrium  and  nonequilibrium  effects. 


mm 


The  discussion  which  follows  has  been  divided  into  four  parts.  Described  in  Section  2. 1  is 
the  steady  inviscid  thermochemical  equilibrium  flow  field  solution  which  has  been  used  to 
provide  the  basis  for  all  of  the  inviscid  plasma  work  done  for  this  report.  Section  2. 2 
describes  the  thermal  equilibrium  multicomponent  chemically  reacting  streamtube  solution 
which  has  been  used  in  conjunction  with  the  equilibrium  flow  field  solution  to  obtain  the 
nmequilibrium  flow  properties.  In  Section  2. 3  the  results,  presented  In  Figures  2.  3-1 
through  2.  3-26  are  discussed.  In  Section  2.  4,  future  flow  field  work  is  discussed. 

2.1  EQUILIBRIUM  INVISCID  FLOW  FIELD  SOLUTION 

The  inviscid  therm ochemical  equilibrium  flow  fields  that  have  been  produced  for  this  study 
were  obtained  from  the  direct  detached  shock  solution  of  Gravalos,  Edelfelt  and  Emmons  1+ 
in  the  transonic  region,  and  by  the  method  of  characteristics  solution  of  Jones  and  Friedhofer2 

3 

as  modified  by  Eds  all  ,  in  the  supersonic  region.  In  the  transonic  region,  the  solution  is 
obtained  in  intrinsic  coordinates  (streamwise  -  normal  coordinates)  by  a  forward-marching 
finite-difference  integration  along  normals.  In  practice,  the  direct  transonic  solution  is 
obtained  by  prescribing  the  body  geometry,  the  free  stream  conditions  and  initial  estimates 
of  the  body  surface  pressure  distribution  and  shock  shape,  solving  the  governing  equations  as 
described  above  and,  then,  based  on  the  results,  modifying  the  surface  pressure  distribution 


♦References  appear  in  Section  5. 


and  shock  shape,  and  iterating  until  the  boundary  conditions  at  the  shock  are  satisfied  to 
within  a  specified  tolerance.  After  this  solution  has  been  determined,  the  results  at  the 
downstream  end  of  the  transonic  region  provide  the  necessary  initial  data  for  the  character- 
istics  solution.  In  both  regions  the  thermochemical  equilibrium  properties  of  air  are 
treated  by  defining  what  is  known  as  a  "y*  gas”  according  to  the  relation: 


in  which  P  signifies  the  static  pressure,  p  the  static  density,  S  the  entropy  and  y*  is  deter- 

4 

mined  by  curve  fits  of  the  data  taken  from  the  NBS  tables  in  the  form: 


y* 


+  b 


(2.1-2) 


where  a  and  b  are  continuous  functions  ol  entropy  and  piecewise  continuous  functions  of 
pressure.  Although  the  y*  gas  is  an  approximate  v  to  real  equilibrium  sir,  its  validity  as 
well  as  the  validity  of  the  methods  of  solution  have  been  substantiated  by  comparisons  with 
other  solutions  and  with  experimental  data. 

2.2  NONEQUIUBRIUM  STREAM  TUBE  SOLUTION 

The  inviscid  thermal  equilibrium,  chemical  non  equilibrium  results  produced  for  this  study 
were  obtained  from  the  streamtube  solution  of  McMenamin  and  O'Brien5.  In  essence,  the 
non  equilibrium  streamtube  method  consists  of  the  finite-difference  simultaneous  solution  of 
the  governing  equations  for  the  steady  inviscid  flow  of  a  reacting  multicomponent  gas  in 
intrinsic  coordinates  by  deleting  the  transverse  momentum  equation(s)  and  independently 
prescribing  the  streamwise  variation  of  one  of  the  flow  properties.  It  should  be  understood 
that  the  streamtube  method  is  approximate  only  to  the  extent  that  the  prescribed  streamwise 
governing  variable  is  approximate,  but  that  to  prescribe  the  streamwise  governing  variable 
exactly  would  require  the  exact  nonequilibrium  multidimensional  fliw  field  solution  to  be 
known.  The  usefulness  of  the  streamtube  method,  therefore,  derives  from  the  fact  +*»»+  for 
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shock  layer  flows  the  static  pressure  is  relatively  insensitive  to  nonequilibrium  effects,  a 
result  which  can  be  ascertained  by  comparing  the  pressurt  fields  for  the  two  limiting  cases 
of  frozen  and  equilibrium  flows.  Consequently,  by  applying  the  nonequilibrium  streamtube 
analysis  to  a  sufficient  number  of  streamlines  from  either  an  equilibrium  or  a  chemically 
frozen  flow  field,  a  good  approximation  of  the  nonequilibrium  flow  field  can  be  obtained. 

Indeed,  when  the  streamtube  method  is  judiciously  applied,  the  results  can  be  exceptionally 
good.  In  detached  shock  layer  flows  the  predominant  nonequilibrium  effect  on  the  shock  is  to 
increase  the  shock  detachment  distance  (by  as  much  as  a  factor  of  2-1/2)  over  the  correspond¬ 
ing  equilibrium  valum  To  a  lesser  degree  the  shock  shape  in  the  nose  region  is  more  blunt, 
lor  the  nonequilibrium  case  than  it  is  for  the  corresponding  equilibrium  case.  When  these  two 
effects  are  taken  into  account  and  me  pressure  that  is  obtained  from  the  thermal  equilibrium 
frozen  shock  flow  solution  is  used  to  replace  the  equilibrium  pressure  behind  the  shock,  the 
exact  boundary  conditions  for  the  nonequilibrium  flow  field  can  be  approximated  quite  closely 
and,  consequently,  the  streamwise  pressure  distribution  used  with  the  streamtube  solution 
can  be  adjusted  even  closer  to  the  exact  nonequilibrium  pressure  distribution  than  if  it  were 
taken  directly  from  the  equilibrium  flow  field. 

As  pointed  out  in  Reference  5,  the  validity  of  the  nonequilibrium  streamtube  solution  has  been 
substantiated  by  comparisons  with  experimental  data  and  with  exact  numerical  flow  field 
solutions  for  the  following  cases: 

a.  In  a  gomparison  with  the  shock  tube  electron  density  measurements  of  Lin,  Neal,  and 
Fyfe  ,  the  results  of  the  steamtube  solution  were  within  the  maximum  estimated 
error  of  the  data  (+50%). 

b.  In  a  comparison  with  the  shock  layer  electron  density  measurements  o'  Kaegi  and 
Chin  ,  the  results  of  the  streamtube  solution  were  within  the  interpretation  of  the 
data  in  the  inviscid  region. 

c.  In  i  comparison  with  the  nonequiiibrium  method  of  characteristics  solution  of 
Gravalos  and  Studerus  for  the  properties  along  a  cone  surface,  the  streamtube 
results  were  in  virtually  exact  agreement. 
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Also,  it  is  of  interest  to  note  that  in  the  electron  density  sensitivity  study  for  sphere-cone 
streamtube  applications,  which  was  reported  in  Reference  5,  it  was  determined  that  for  a 
constant  difference  of  +10%  in  the  streamwise  pressure  distribution  (which  is  representative 
of  the  largest  error  introduced  by  assuming  an  equilibrium  pressure  distribution)  the 
electron  density  was  affected  by  less  than  +50%. 

2.3  RESULTS 

The  configuration  selected  for  this  study  is  an  8. 6  degree  sphere-cone  with  a  nose  radius  of 
0. 5  inch.  According  to  the  method  of  solution  discussed  in  Section  2. 1,  flow  fields  for  this 
configuration  were  obtained  at  altitudes  of  50,  100,  150,  200  and  250  kilofeet.  The  free 
stream  conditions  for  each  of  these  flow  fields  are  listed  in  Table  2. 3-1.  *  The  first  four  of 
these  flow  fields  were  determined  under  the  condition  of  chemical  equilibrium,  whereas  at 
250  kilofeet  a  perfect  gas  (y  =  7/5)  flow  field  was  obtained  to  serve  the  double  purpose  of 
providing  steamwise  pressure  distributions  needed  in  the  nonequilibrium  streamtube  analysis 
for  the  plasma  attenuation  study,  and  als<  •  to  provide  the  gas  dynamic  properties  in  the 
limiting  case  of  no  vibrational  excitation  which  will  be  needed  in  the  antenna  breakdown 
analysis. 

The  reaction  system  used  in  the  nonequilibrium  streamtube  solution  along  with  the  kinetic 
rate  coefficients  and  the  collision  efficiencies  is  shown  in  Table  2.  3-2.  The  nonequilibrium 
results  which  are  presented  in  Figures  2.  3-1  through  2,  3-26  were  obtained  in  the  following 
manner.  Using  streamwise  pressure  distributions  from  the  perfect  gas  or  equilibrium  flow 
fields  for  selected  streamlines  (which  herein  are  identified  by  the  bow  shock  angle  at  the 
streamline-bow  shock  intersection),  the  nonequilibrium  streamtube  solution,  wnich  is  coded 
in  Fortran  IV  and  run  on  the  IBM  7094,  is  obtained  for  one  streamline  at  a  time,  and  the  flow 

*In  the  interest  of  brevity,  floating  point  numbers  in  Table  2. 3-1  through  2. 3-7  have  been 
expressed  by  means  of  a  slash  (after  and  slightly  above  the  last  significant  figure)  followed 
by  the  pow  3r  of  10.  For  example,  the  number  1. 234  x  10-8  is  written  1. 234  V8.  in 
addition  it  is  to  be  understood  that  for  the  columned  numbers  the  came  exp  onent  also  applies 
to  any  number  below  written  without  one. 
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Table  2.3-2.  Reaction  System 


Number 

1 

°2  +  M 

0  +  0  +  M 

2 

N2  +  M 

N  +  N  +  M 

3 

NO  +  M 

N  +  0  +  M 

4 

NO  +  0 

02  +  N 

5 

n2  +  0 

NO  +  N 

6 

N  +  0 

NO++  e' 

>v.  • 


* 


00,N  ,  N,  O,  NO 


Reaction  Rate  Consiants,  k 


B„ 


exp  -(C,  /RT) 


2.5\16 

2.0\21 


8. 9\14 
1.9l\20 


5.5\20 

3.2\9 


1.67\20 

2.67^° 


0.8\13 

6.4V 


2.X13 

7.0\18 


B 

r 

cf 

■Cr 

-.44 

117945 

0 

-1.57 

225014 

0 

* -1.52 

150002 

0 

.92 

39100 

7042.5 

-.04 

75011.5 

0 

-.9 

64356.5 

0 

NOTE:  Units  are  moles/cc,  sec. ,  K,  and  cal/gmol 
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properties  are  stored  at  the  intersections  of  the  streamline  with  the  specified  body  normals. 
After  the  last  streamtube  run  is  made,  the  position  of  each  streamline  along  each  normal  is 
determined  by  continuity,  using  the  nonequilibrium  values  for  the  density  and  velocity.  By 
selecting  a  sufficient  number  of  streamlines  and  body  normals,  the  end  result  is  a  well 
defined  nonequilibrium  flow  field  which  is  approximate  only  to  the  extent  to  which  the  bow 
shock  and  the  streamwise  pressure  distributions  are  approximate. 

For  the  present  results,  no  attempt  was  made  to  adjust  the  perfect  gas  or  equilibrium  shock 
shapes  for  nonequilibrium  effects.  However,  the  pressure  immediately  behind  the  bow  shock 
was  in  every  case  replaced  by  the  value  obtained  from  the  thermal  equilibrium  frozen  shock 
flow  solution.  Although  the  shock  shape  was  not  adjusted,  it  is  significant  to  note  that  after 
constructing  the  approximate  nonequilibrium  flow  fields,  a  comparison  showed  that  the 
maximum  percent  difference  in  the  shock  radius  is  -3.  5%  between  the  nonequilibrium  and 
frozen  values  and  +3%  between  the  nonequilibrium  and  equilibrium  values.  In  Table  2.  3-3, 
along  with  a  summary  of  the  body  geometry  corresponding  to  three  normals  which  have  been 
selected  to  display  the  nose  bluntness  and  altitude  effects  on  the  nonequilibrium  flow 
properties,  a  comparison  is  presented  of  the  shock  radii  at  the  intersections  of  Ihe  shock 
and  the  normals  as  determined  by  the  different  solutions.  By  far,  the  larger  part  of  the 
error  is  due  to  the  fact  that  the  perfect  gas  and  equilibrium  flow  field  values  for  the  shock 
detachment  distance  were  not  adjusted  for  nonequilibrium  before  being  used  with  the  stream- 
tube  solution.  This,  however,  will  not  significantly  affect  the  results  even  along  the  normal 
on  the  noue. 

The  altitude  trends  of  the  shock  radius  comparison  in  Table  2.  3-3,  particularly  for  the  first 
normal,  are  significant  and  deserve  further  attention.  It  may  be  noted  that  at  50  kilofeet  the 
values  of  the  shock  radius  according  to  the  equilibrium  and  nonequilibrium  solutions  agree  to 
within  0. 06%,  which  indicates  that  the  flow  is  nearly  in  equilibrium.  As  the  altitude  increases, 
the  nonequilibrium  value  for  the  shock  radius  becomes  increasingly  larger  than  the 
equilibrium  value,  indicating  that  the  flow  is  becoming  further  out  of  equilibrium.  The 
comparison  at  250  kHofee1:  with  the  perfect  gas  value,  however,  shows  that  even  at  this  altitude 
the  flow  is  not  frozen  in  the  nose  region. 


Table  2.3-3.  Summary  of  Geometric  Variables  for  Selected  Body  Surface  Normals 
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In  Figures  2.  3-1  through  2.  3-10  the  contrast  between  chemical  equilibrium  and  non¬ 
equilibrium  effocts  at  an  altitude  of  150  kilofeet  is  demonstrated  by  plotting  the  streamwise 
variation  of  the  flow  properties  against  the  distance  along  the  body  surface  for  three  stream¬ 
lines  which  intersect  the  bow  shock  at  angles  of  approximately  50,  70  and  89  degrees.  In 
sequence,  these  plots  show  this  variation  for  the  pressure,  temperature,  density,  velocity 
component  parallel  to  the  surface  (all  normalized  by  the  corresponding  free  stream  quantities), 
and  the  number  densities  of  molecular  oxygen,  molecular  nitrogen,  atomic  nitrogen,  atomic 
oxygen,  nitric  oxide  and  electrons.  Before  proceeding  to  a  brief  discussion  of  some  of  the 
more  interesting  aspects  of  these  results,  it  should  be  pointed  out  that,  in  Figure  2.3-1,  due 
to  the  positive  normal  pressure  gradient  in  the  nose  region,  the  higher  shock  angle  stream¬ 
lines  initially  expand  more  rapidly;  however,  in  the  over-expansion  region  this  trend  is 
reversed  prior  to  converging  at  the  cone  surface  pressure  value. 

For  the  nonequilibrium  results  shown  in  Figures  2.  3-1  through  2.  3-10,  it  may  be  noted  that 
the  flow  on  the  89  degree  streamline  is  almost  in  equilibrium  at  the  first  normal.  This  is 
caused  by  its  passage  through  the  near  stagnation  region.  At  this  point  oxygen  is  almost 
entirely  dissociated.  As  the  flow  continues  to  expand  over  the  nose,  atomic  oxygen  freezes 
in  the  nonequilibrium  case,  whereas  for  equilibrium  it  begins  to  recombine.  Ibis  single 
feature,  more  than  any  other,  determines  the  nonequilibrium  effects  which  are  indicated  in 
these  results. 

In  Figures  2.  3-11  through  2.  3-25,  nose  bluntness  and  altitude  effects  on  the  nonequilibrium 

flow  properties  are  shown  by  plotting  the  flow  properties  against  shock  layer  thickness,  with 

altitude  as  a  parameter,  for  three  selected  normals  located  along  the  body  surface  at  s,  /R.t 

b  N 

equal  to  0. 4962,  61. 25  and  243. 3.  Sequentially  in  groups  of  three,  the  normal  profiles  are 
given  for  the  pressure,  temperature,  density  and  velocity  parallel  to  the  surface  (all 
normalized  by  the  corresponding  free  stream  quantity)  and  the  electron  number  density. 

Profiles  of  the  remaining  species  number  densities  are  presented  in  Table  2.  3-4,  2. 3-5  and 
2. 3-6  for  the  different  altitudes  at  the  selected  normals.  In  Figure  2. 3-26  the  integral  of  the 
electron  density  across  the  shock  layer  is  plotted  against  the  normalized  coordinate  along 
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NORMALIZED  VELOCITY  PARALLEL  TO  SURFACE 


NORMALIZED  COORDINATE  ALONG  BODY  SURFACE,  8,/*^ 


Figure  2. 3-4.  Streamwise  Variation  of  Normalized  Velocity  Parallel  to  Surface  vs 

formalized  Coordinate  Along  Body  Surface 


2-13 


•*;4Kia*ki«a 

isntiitt 

&f4HSABI| 


_ 1 

‘.'ii'jrz  L  jtt 

T 

■ 

L_5 

"  150  KUO 

CM 

war 

'•T. 

Zgq 

1  " 

■MOCK  A 
X)  (DEC 

NGLK 

2) 

:::::::::::::::::::: 


HiSIls 


iui  IISZ2XS*  :  = 

::i!H{i!!!ll— .  t  SS3H  SSii 


!j!ini>pifiii!ftil  iitljiifiiiiij 


S&iiififflifii 


w.iKMinfin 

)§'  »M> 


«n  iiMniimn  <•••  •••«! 

•u  !«•■ 

hi  until*  iM(:illN9fla( 

-  iii* '<ifi  iiiii  liar 
*i>  .<n ''.iiiiiii iiiii 

•ii 

.1 ;  >1111  mil  mil  ■%««  III*;: 

ll>  Mll»H;:-:(l|||iMj|||| 


IllUlSt 


iMiiiiia 

■  ■■ 

■ 

■i 

:::::::::::::: 

till . . 

: 

■ 

■ 

Mil  Nil 

UN 

* 

■ 

MU  Nil 

!■■■■ 

a 

(do/'oo)  nn  'MaooH.tiN  HvinoaaoH  jo  ausn3q  aaawnN 


IftHSlEl&ffMiM! 


Figure  2. 3-7.  Streamwise  Variation  of  Atomic  Nitrogen  Number  Density  vs  Normalized 

Coordinate  Along  Body  Surface 
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Igure  2. 3-10.  Streamwise  Variation  of  Electron  Number  Density  vs  Normalized  Coordinate 

Along  Body  Surface 
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Figure  2. 3-15.  Normalized  Temperature  vs  Shock  Layer  Thickness;  sfe/RN  =  6. 125\1;  Nonequilibrium 
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Figure  2. 3-16.  Normalized  Temperature  vs  Stock  Layer  Thickness;  8,/Ii  =  2. 433 \  ; 

Nonequilibrium 
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ELECTRON  NUMBER  DENSITY,  N  (no./cc) 


0  0. 02  0. 04  0. 06  0. 08  0. 10 


SHOCK  LAYER  THICKNESS,  YAg^ 


Figure  2. 3-24.  Electron  Number  Density  vs  Shock  Layer  Thickness;  s^/RN 


6. 125^ 
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Table  2. 3-4.  Altitude  Dependence  of  the  Species  Number  Density  Normal  Profiles 

S.  /R  =4.962  “1 


Species  Number  Density  (particles  per  cc) 


Altitude 

Y/Y 

o9 

No 

(kilofeet) 

sn 

2 

2 

50 

2. 723X1 

5.701^° 

6.520 

2. 2. 8^\9 
2.37! 

I:Joo\° 

9.689 

6. 662 

2.651 

2.509 

100 

7.  74  V 
2.429V 
4.838 

8.  29oV 
8.742 
9.556 
1.374V 

2.  05X8 

2.  080 

2. 1584 

2.  279 

7.005  0 

1. 000  \ 

3.100 

6.399V 

2. 384 
2.410 

150 

8.06X1 

7.4oX3 

1.12X4 

2.  097V 
2.000 

2.393^ 

2.655 

2. 127 

4. 780 

7.244 

1.910^ 

6.717\16 

2.187 

7,044  0 
1.00  X 

2.238 

2. 530 

200 

7. 98qV 

5.31X8 

1. 226V 

3.  072s}6 
3.040 

2.706N1 

3.068 

4.782 

3.220 

7.  098 

1.  ooox 

l.‘027\16 

3.463 

3.867 

250 

4.89^ 

1.975N 

1.  519s}4 
3.244 

6.774 

3. 73lV 
3.845 

4.101 

4.612 

9. 954  15 

1. 192\ 

4.412 

7.594 

4.693 

1.000 

1.384 

5.211 

6.78^s8 

6.405 

5.228 

0 

1.17lX8 

1.707 

1.660 

1.503 

1.195 

0 

2.  Ol^J7 

1.838 
1.484 
1. 103  i_ 
7.386^® 

0 

2. 194s\6 
1.707 
1.145  ie 
6.54^5 
2.973 
0 

8.  299^4 

5.944 

3.517 

U  80ll3 

6.52RX3 

0 


1. 296s\9 
1.329 
1.418 
0 

1. 51^\8 

1.526 

1.553 

1.566 

1.510 

0 

1.617\7 
1.562 
1.452 
1.313 
1. 113 
0 

2. 080^  6 
1.828 
1.541 
1.324 

8.  OO2V 

0 

1. 522s^5 
1.252 
7. 516s?4 
3.348 
1.062 
0 


1.12lV 

1.203 

1.476 

0 

3.8i(V 
3.867 
4.161 
5. 2P3 
8. 962 
0 

3. 753s\5 
4.945 
8.684 
1.693V 
3.174 
0 

1.624s}5 

2. 878 

4.567 

4.656 

2.588 

0 

5.1oH4 

4.070 

2.238 

8.84<Kr 

2.459 

0 


"W- 

■  I  . 


Table  2. 3-5.  Altitude  Dependence  of  the  Species  Number  Density  Normal  Profiles, 

s,/R  =  6.125  V 
b  N 


(ldlofeet) 


Species  Number  Density  (particles  per  cc) 


0  '-3 

4.  213'V 

1.247\2 


2.416 
4.  088 
5.726 
7.875 


1.617' 
2. 259 
5.303 
1.366 
2. 873 
5.636 


1.647H8 


1  •  u  1  v  « 

1. 145V 

1.  000\ 


1.  094 
1.847 
3.413 


1.69S 

1.  844 

2.  018 
2.469 
3. 169 
4.138 
6.955 


2.  962s 
2.  214 


9.  882\ 
3.!02  14 


L.285V 


8.65l\n 

1.614N1" 

2.339\ 


7. 956s\7 
8. 009 
7.917 
6.872 
6.238 


4.946^" 
1. 39J^\ 

0 


0 

3.832 

8. 272  _ 

1.279^ 

1.904 

2.439 

4.013 

5. 907 

8.275  . 

1.166V 

2. 214 

1.000 


l.ioiM2 

1.206 

1.351 

1.706 

2.582 

4.354 

5. 815 

4.500vl6 

9.934 


L.  57^\7 


v  I  «7VT  -  — 

1.615V 


3.493 

3.330 


1.636 
1.739 
1.853 
2.048 
2.  231 
2.558 
2.898 
3.763 
6.082  ,c 
1.315^8 
1.254 


8. 295  \16 
8.057 
7.733 
6.  865 
5.433 


3.844 


3. 992' 


1.062 

1.082 

1.128 

1-166 

1.231 

1.286 

1.184 


4-^o 

2  417' 
5.249\2  . 
1.7518\15 
0 


1  A 

4.184X 

1.733V 
1. 005V 


0  4 
7.858V 
2. 501V 


5. 777 
9.369  „ 
1.484V 
2.  002 
3.479 
5.117 
7. 224 
1.  03(V 
l.OO^s! 


2.516^° 

2.190 

1.809 

2.169 


5.688  „ 
1.64^5 


6.374 
1. 032N\6 
1.674 
3.381 


1.  99X6 

1.931 

1.912 

1.924 

1.945 

1.984 

2.008 

2.262 

2.721 

3.890 

6.306  ... 
_ _ vi7 


1.  843s \6 


L.  550^ 


1.624 

1.270 


1.273X 


9. 21  <K 
6.289 
2.900  .  . 

7'16Ci3 

4  I9i>  ** 

i:«^2 

5. 873)! 
4.82l\° 

0 


1. 55 0\ 

1.457 

1.352 

1.266 

1.199 

1.128 

1.057 

6.532^" 

1>14ll3 

1.324V 


1.793^ 


7.437s\6 

8‘42V7 

1.144V 

1.523 

1.535 

2.473 

5.603V5 

1.362X 

0 


7.6R0\ 

7.306 

6.396 

6.684 

7.097 

9.189 

3.00lV 

6.632xM 

1.936\‘0 

1.49^ 

0 


4.34^\ 
4.315 
7. 737  , 


2. 341V\2 


7.363 


s.wsV 


3. 083V 
2.999V" 

7*51^13 


1.199' 


2. 006X 


«  as 


’  * 

V  V  -t  v.« 


■ 


'*■  ?: 


Table  2.3-5.  Altitude  Dependence  of  the  Species  Number  Density  Normal  Profiles, 

su/R.=6.125\1  (Cont) 

D  XM 


Altitude  Y/Y 
(Idlofeet) 


Species  Number  Density  (particles  per  cc, 


.  : 


9.506\1 
3. 69<N"3 
7.235  ’ 
1. 17sV~ 
1.821 
2.381 
3. 932 
5.  m 
7.596  1 

1.  000 X 


1.  379^® 
2.14lV 

4.829\ 

l.*5269\13 
8.078  .. . 
2.284\14 
7. 121 
1. 126\15 
1.568 
2.695 
5.  256 


2.  789^5 
2.655 
2„  bo9 
2. 543 
2,  557 
2.674 
2.809 
3.288 
4. 305 
5.908  .. 
!.  015\16 
1.980 


i.78(^5 

UV 

6.  763\ 
3.753 
2.199 
7. 925N\3 
2.821  10 
1.830\3 
2.  864\J° 

5.450V 
2. 168\ 

0 


1.953^5 
1.728 
1. 523 
1.393 
1.239 

8.  013V 
2.365 

2.3m.: 

4.  08$V 
3.47l\6 
0 


2. 712^2 

iSfc* 

6.783 

1.982\14 

3.214 

3.279 

1.164 

1.435\J? 

3.617V1 
3. 877 v 
0 


*  * 


3. 088V4 
3.543VJ 
9.200 

1.668V 

2.511 
3.639 
4.668 
7. 038 

9.495  , 
1.411X1 

2.495 
l.OO^V 


2.575 

4.199 

5.486 

6.556 

7.713 

9.840 

1.286\W 

1.631 

2.745 

6.962 

7.302 


2. 712^4  suit's3  1.223V4 


2.719 

2,711 

2. 76  U 

2.869 

3.010 

3.234 

3.971 

4.902 

6.149 

1.034\15 

2.622 

2. 750 


3.993 
2.  228 
1.209  19 
6. 996\ 
4.101 
2.060  „ 
9.455V11 
9.  069^ 
2.478\® 

3. 099\ 
6.678V15 


1.068 

7.805V13 

5.188 

3.460 

3. 309 

1.394 

1.086 

2.382V12 

2. 902V1 

1.696V 

1.179V10 


3.509^3 

3.411 

2.764 

1.959 

^112 

9.12lV 

5.583 

4.720 

1.163  ,, 

1.852V11 

1. 805X 

2.969V 

0 


•’4* 
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Table  2, 3-6.  Altitude  Dependence  of  the  Species  Number  Density  Normal  Profiles 
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the  body  surface  with  altitude  as  a  parameter.  Some  of  the  more  noteworthy  aspects  of  these 
results  are  discussed  in  the  following  paragraphs. 


It  is  worth  noting  in  Figures  2.  3-11,  2.  3-12  and  2.  3-13  that  the  pressure  profiles  obtained 
i  om  the  perfect  gas  or  equilibrium  flow  fields  actually  follow  the  free  stream  Mach  number 
variation  rather  than  the  altitude  variation.  Also,  for  these  pressure  profiles,  which  are 
derived  from  the  perfect  gas  or  equilibrium  streamwise  pressure  distributions  and  reproduced 
during  the  construction  of  the  approximate  nonequilibrium  flow  field,  the  pressure  behind  the 
shock  is  replaced  by  the  value  from  the  thermal  equilibrium  frozen  shock  flow  solution.  This 
change  at  Y/Yg^  =  1  is  particularly  noticeable  on  the  first  normal,  where  the  result  is  a 
lower  value  for  the  equilibrium  cases  and  a  higher  value  for  the  perfect  gas  case. 

Some  attention  should  also  be  directed  to  the  development  of  an  off-wall  peak  in  the 
temperature  profiles  back  on  the  conical  surface  (Figures  2.  3-15  and  2.  3-16)  as  the  altitude 
is  increased.  This  trend  is  unquestionably  one  of  the  more  marked  effects  of  nonequilibrium 
flow,  inasmuch  as  the  temperature  for  equilibrium  flow  decreases  monotonically  from  the 
wall  to  the  shock  (except  in  the  region  where  the  streamlines  have  traversed  that  portion  of 
the  shock  where  the  minimum  shock  angle  occurs  as  a  result  of  the  overexpansion). 

2.4  FUTURE  FLOW  FIELD  WORK 

As  required  for  use  in  the  antenna  breakdown  study,  inviscid  flow  properties  at  several 
altitudes  above  175  kilofeet  will  be  obtained  for  the  limiting  cases  of  no  vibrational  excitation 
and  full  vibrational  excitation,  for  critical  comparison. 

Chemically  reacting  boundary  layer  calculations  for  ablation  product  contaminated  flow, 
including  the  important  metallic  trace  impurities,  will  be  made  at  various  altitudes  to 
determine  viscous  and  contamination  effects  on  plasma  properties. 

Furthermore,  a  survey  of  the  current  state  of  the  art  in  body  flow  field  analysis  will  be  made 
and  reported  late  in  the  contract  period. 
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SECTION  3 

ANTENNAS  AND  PROPAGATION 

\ 

3.1  ANTENNAS 

3.1.1  ECM  ANTENNAS 

This  subject  is  covered  in  entirety  in  Volume  II  of  this  report. 

3.1.2  ANTENNA  FIELD  THEORY 

The  purpose  of  this  discussion  is  to  survey  the  methods  of  calculating  antenna  near-field 
distributions,  which  should  be  known  before  antenna  breakdown  and  reentry-induced  plasma 
effects  can  be  analyzed. 

In  general,  the  complete  radiated  fields  of  a  given  antenna  must  be  calculated  in  terms  of  the 
antenna  geometry,  current  distribution,  and  wavelength.  For  prediction  of  the  effects  of 
reentry  plasmas  on  the  antenna  properties  and  patterns,  it  will  be  necessary  to  consider  the 
antenna  as  a  finite,  continuous  source.  Thus  calculation  of  the  transmitted  power,  patterns, 
and  polarization  becomes  a  solution  of  an  exact  boundary-value  problem.  Htwever,  since  the 
wave  equation 

Vx  VxE  =  p  6  (r  -  a) 

is  separable  in  only  a  few  restricted  coordinate  systems  (e.  g. ,  spheroidal  and  infinite 
cylindrical),  approximate  methods  usually  must  be  employed. 

A  common  method  is  an  integral  equation  approach.  From  Huygen's  principle  we  can 
evaluate  the  fields  external  to  a  surface  by  the  specification  of  the  charges  and  currents  on 

"4 

that  surface.  In  terms  of  the  magnetic  vector  potential  A, 


VxA, 


(3. 1.2-1) 
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E  =  -  V  x  ^  x  A 

For  a  line  source  of  length  L, 


^o 

4  IT 


(3. 1.2-2) 


(3. 1. 2-3) 


where  the  integration  is  over  the  source  coordinates  (xf,  y',  z’)  and  we  have  assumed  the 
source  current  I(z')  along  the  zf  axis.  Now,  in  principle,  the  Poynting  vector,  S  =  ExH,  is 
the  solution  of  our  problem,  specified  by  Equations  3. 1. 2-1  through  3. 1. 2-3.  However,  the 
integral  (Equation  3. 1. 2-3)  must  be  evaluated  approximately.  It  is  convenient  to  consider  the 
fields  in  three  different  regions; 


r  < 


Fresnel  zone, 


2L 

r  >  —  ,  far  or  radiation  zone. 

If  the  source  were  a  surface  or  a  volume,  then  the  Integral  in  (Equation  3. 1. 2-3)  would  be 

modified  accordingly.  The  essential  problem  is  the  determination  of  the  current  distribution, 

I  (z').  The  general  theory  in  terms  of  arbitrary  source  current  distributions  has  be6n  given 
9 

by  Silver  .  We  will  be  mainly  interested  in  the  solutions  for  dipoles  (finite  cylindrical 
antennas)  and  for  traveling-wave  antennas  (helices).  The  procedures  required  in  the  analysis  of 
these  two  types  can  then  be  applied  to  other  types  of  antennas. 
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Traveling-wave  antennas  may  be  represented  by  letting  the  current  in  Equation  3. 1. 2-3 
become  sinusoidal,  so  that 

OO 

I  (z»)  =  E  I  exp  C-j  {n7rz’/L}]  (3. 1.2-4) 

n=oo 

z  +  L 

I  (z*)  exp  [j  (nTrz’/L)]  dz'  (3. 1. 2-5) 

Thus  the  source  is  represented  as  a  sum  of  traveling  waves.  The  near  fields  of  a  line  source10 
represented  in  this  manner  cau  be  e  valuated  by  quantizing  th..  length  into  segments  so  short 
that  the  field  point  is  in  the  far  field  of  the  segment  while  still  in  the  near  field  of  the  whole 
antenna.  Thus  the  integral  of  Equation  3. 1. 2-3  is  broken  up  into  N  integrals,  in  each  of 
which  the  current  I  is  a  constant.  The  integration  is  simplified  ^othat  the  near-field  is  of  the 
form 
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The  same  method  may  be  generalized  to  rectangular  and  circular  sources. 


An  important  type  of  traveling-wave  antenna  is  the  helix,  which  is  a  slow-wave  structure. 
Using  the  sheath  helix  approximation,  this  problem  can  be  solved  exactly11’12. 


A  complete  discussion  of  the  properties  of  single-wire  and  multiwire  helices  is  given  by 

Watkins  with  impedance  calculations,  k  —  /9  diagrams,  and  field  patterns.  In  addition, 

Walter10  discusses  helices  in  the  context  of  traveling-wave  structures.  The  effects  of 

ionized  atmospheres  on  the  transmission  line  characteristics  of  helices  has  been  discussed 
14 

by  Neurather,  et  al. 


The  dipole  or  cylindrical  antenna  has  been  analyzed  thoroughly  by  King  .  In  order  to 

calculate  the  radiation  patterns,  impedances,  and  near  fields  of  finite  cylindrical  antennas, 

is 

it  is  necessary  to  solve  the  well-known  Hallen  integral  equation  .  This  was  done  using  an 

17 

iterative  method  assuming  a  sinusoidal  current  distribution,  as  outlined  in  Kraus  .  We  then 

obtain  well-used  formulas  for  the  impedances  and  near-fields  as  a  function  of  wavelength  and 

18 

antenna  dimensions.  A  more  accurate  determination  has  been  made  by  King  and  Middleton 

assuming  a(sinO)/®  initial  current  distribution.  Over  the  years,  a  voluminous  literature  has 

19  20  21 

been  generated  on  the  subject  '  ’  .  The  essential  problem  is  determining  the  axial 

current  distribution;  the  calculation  of  the  fields  and  impedance  follows  directly  from  the 
methods  given  in  References  15  and  17.  Since  the  Initial  current  distribution  is  so  important 

there  can  be  no  true  separation  from  plasma  effects.  The  analysis  of  dipoles  in  ionized 

22 

atmospheres  has  been  performed  by  Balmain  ,  who  assumes  a  triangular  current  distribution, 

23 

and  by  Cook  and  Edgar  ,  who  solve  for  the  current  distribution  assuming  that  the  electron 

velocity  vanishes  at  the  edge  of  the  plasma.  For  the  frequency  and  power  ranges  under 

consideration,  dipole  (or  multipole)  radiators  are  more  important;  however,  calculations  for 

24 

loops  similar  in  method  have  been  given  by  Iizuka  and  King  . 


3.2  LINEAR  PLASMA  EFFECTS 

As  indicated  in  Section  3. 1. 2  on  antenna  field  theory,  the  near  fields  and  radiated  fields  of  a 
radiating  source  are  determined  by  the  current  and  charge  distributions  of  the  source.  The 
presence  of  a  plasma  layer  will  affect  the  current  distribution  and  impedance  of  an  antenna 
since,  basically,  the  plasma  dielectric  constant  differs  from  that  of  free  space, 


(o  (w  *"  Jv  c)  ’ 


(3. 2-1) 
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(3. 2-2) 


where 


v  =  Electron  collision  frequency 

n  =  Electron  density 

e  =  Electron  charge 

m  =  Electron  mass 

Thus  the  wave  number 

27T  2tt  V  1 

K  ~  \  ~  c  ~  c  ’  °  =  r  (3-2-3) 

Jn  € 

woo 

and  propagation  constant,  y,  must  be  calculated  in  order  to  determine  the  field  configurations. 
However,  the  presence  of  the  plasma  layer  allows  the  existence  of  many  modes  of  propagation 
in  addition  to  the  continuous  spectrum  of  free-space  wavelengths.  Therefore,  the  analysis  is 
made  convenient  by  first  considering  the  possible  modes  in  a  simple  structure,  e.  g. ,  a  plasma 
layer  over  a  plane  surface,  and  then  considering  other,  more  practical  geometries. 

The  total  field  of  this  antenna  configuration  can  be  related  to  the  impedance  and  dimensions  of 
the  source  antenna  and  the  plasma  layer  by  integral  transform  techniques.  Applying  the 

Huygens  theory  of  Equation  3. 1. 2-1,  3. 1. 2-2  and  3. 1. 2-3. ,  we  are  led  to  the  evaluation  of  the 
integral, 
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+r 
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exp  (-yz)  dz 
(3.2-4) 


where  the  propagation  constants  in  the  z-direction, 


2  2  1/2 
Kx  =  (k2  e  +  y)  , 


(3. 2-5) 


2  2  1/2 
Ko  =  (k  +  y)  , 


(3. 2-6) 


and  the  impedance  of  the  plasma  layer  is 


r  = 


Kq  +  )  cot  Ki  d 

K  -  j  K  cot  K  d 
o  1  1 


(3.2-7) 


Thusy  evaluation  of  the  field  E  depends  on  the  location  of  the  singularities  in  the  integrand  of 

y 

Equation  3. 2-4.  If  we  make  the  transformation 


y  =  a+j/3  -*  j  k  p  =  j  k  (£  +  j  r\) 


(3.2-8) 


the  steepest  descent  contour  C  of  Equation  3. 2-4  is  shown  in  Figure  3. 2-1.  The  location  of 
the  poles,  given  by  Equations  3. 2-5,  3. 2-6,  and  3. 2-7,  relate  the  plasma  properties  to  the 
antenna  fields.  The  formulation  may  again  be  generalized  to  two-  and  three-dimensions. 


The  effect  of  a  plasma  sheath  on  an  antenna  has  been  studied  by  many  authors.  Wait  has 
analyzed  the  effect  of  a  thin  plasma  sheet  on  a  dipole  antenna.  He  uses  a  direct  solution  of 
the  boundary-value  problem  in  order  to  calculate  the  radiation  patterns.  He  considers 
electron  collisions  and  obtains  results  with  and  without  an  impressed  magnetic  field.  He  nat 
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Figure  3. 2-1.  Steepest  Descent  Contour  C 


formulated  the  problem  in  terir  s  of  the  radial  Hertz  potential  and  thus  has  obtained  a 
splitting  of  the  field  into  a  radiation  field  (saddle-point  contribution)  and  surface  field  (pole 


contribution).  The  same  configuration  has  been  treated  in  great  detail  by  Tamir  and  Oliner  , 

27  28  10 
Wenger  ,  Bernard  and  Ishimaru  ,  and  Walter  . 


The  propagation  constants,  K,  and  impedances,  T,  may  be  calculated  by  variational  methods  , 
Thus  the  change  in  propagation  constant,  8 ,  with  variation  in  plasma  dielectric  constant,  e,  is 
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The  variational  method  can  also  be  used  to  determine  the  effects  of  slots  on  the  propagation 
constant  of  a  waveguide  or  radiating  structure,  which  is  useful  in  determining  the  near-field 
patterns  for  the  previous  section. 

References  28  and  10  present  the  effects  on  the  radiation  patterns  of  the  various  pole 
contributions.  The  classification  into  pole  contributions  (leaky-wave,  backward  and  forward 
waves,  surface  waves,  etc. )  is  useful  An  studying  the  myriad  of  possible  propagation  modes 
discussed  in  References  26  and  27. 

The  calculation  of  the  dipole  fields  in  the  presence  of  a  plasma  layer  has  been  done  very 
completely  by  Banos.  He  gives  complete  formulas  and  radiation  patterns  for  the  case  of 
an  infinitesimal  (Hertz)  dipole. 

An  alternate  representation  is  obtained  by  angular  transformation  techniques  rather  than 

normal  mode  formulations  based  upon  orthogonal  coordinates.  Thus  the  residue  series 

31  32 

obtained  by  a  Watson  transformation  ’  can  be  evaluated.  This  has  been  discussed  by 

33  34 

Wait  and  by  Bremnvjr. 

The  linear  effects  of  plasma  sheaths  have  been  calculated  for  the  spherical  case  by  Wait  and 
35  36 

Spies  and  for  the  spheroidal  case  by  Yeh.  These  two  particular  geometries  can  be 

solved  by  exact  boundary- value  methods  and  thus  provide  convenient,  although  approximate, 

calculations  of  linear  plasma  effects.  The  slotted  cylindrical  antenna,  which  is  a  useful 

approximation  to  the  reentry  vehicle,  has  been  studied  in  great  detail  by  many  authors, 

37  38  39 

particularly  Smith  and  Golden,  Swift,  and  Rusch.  The  cylindrical  geometry  is 

important  when  the  wavelength  is  comparable  to  the  vehicle  size,  for  then  the  curvature  of  the 

vehicle  and  the  surrounding  antenna  must  be  considered.  Since  the  wave  equation  is  separable 

in  cylindrical  coordinates,  this  is  a  suitable  mathematical  model  which  simulates  the  aft 

portion  of  many  reentry  vehicles.  The  finite  slot  should  be  the  appropriate  representation  of 

the  source,  but  this  involves  complicated  mathematical  procedures  which  restrict  complete 

pattern  calculations  to  plasma  sheaths  that  are  homogeneous  and  loss-less8  ,40*41.  The 

42 

problem  of  the  homogeneous,  lossy  plasma  has  been  solved  by  Sengupta  who  included 
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collisions  through  the  use  of  a  thin  sheath  and  high  and  low  frequency  approximations. 

However,  the  practical  solutions  of  these  equations,  in  general,  are  not  specifiable  because 

the  plasma  properties  may  vary  arbitrarily  within  the  shock  layer.  Therefore,  analytical 

techniques  or  WKB  approximations  are  not  realistic  for  many  problems  of  interest  and  must 

be  abandoned  in  favor  of  exact  numerical  techniques.  The  inhomogeneous  case  is  solved  by 
43 

Tyras  and  Harley,  who  consider  the  plasma  as  a  multilayered  series  of  lossy,  homogeneous 
plasma  slabs. 

This  approach  requires  application  of  the  boundary  conditions  a  great  many  times  and  may 

lead  to  computational  difficulties.  An  alternate  scheme  is  to  integrate  the  wave  equation 
44  45  38 

directly.  ’  Swift  explains  this  method  in  detail  and  gives  antenna  patterns  for 
homogeneous  sheaths  of  arbitrary  thicknesses  and  with  vjw  varying  from  0  to  0.  30. 

The  effects  of  a  conical  plasma  sheath  are  important  for  nose  antennas  where  the  boundary 

layer  is  ionized.  Jince  the  wave  equation  is  not  separable  in  these  coordinates,  various 

mathematical  devices  have  been  applied.  Yeh  uses  an  application  of  Sommerfield's  complex 

order  wave  functions  to  obtain  the  far-field  radiation  patterns.  His  analysis  is  especially 

applicable  to  high  frequency  (ka  >  3)  cases.  In  addition,  his  formulation  of  the  problem  allows 

consideration  of  an  infinite  sheath  whose  thickness  and  electron  density  vary  linearly  with 
w  47 

radial  distance.  Banos  et  al  have  given  a  complete  analysis  of  the  radiation  pattern 

distortion  by  an  infinite  conical  plasma  sheath  around  an  electric  dipole.  Using  the  Watson 

transformation,  they  obtain  a  residue  series  representation  of  the  radiated  fields.  This 

analysis  is  useful  since  it  separates  the  fields  into  surface-wave  and  dipole  contributions. 

The  plasma  model  used  was  collisionless  and  homogeneous.  Jordan  et  al  have  extended  the 

49 

analysis  to  include  electron  collisions.  •  Pridmore  Brown  has  used  the  Wiener  Hopf  technique 
to  extend  the  Banos  solution  to  the  finite  length  cone.  Again  these  analyses  hold  for 
homogeneous  thin  plasmas. 

Finally,  we  may  note  the  utility  of  the  linear  plasma  solutions  in  obtaining  nonlinear  calculations. 
50 

Menzel  gives  a  general  discussion  of  the  construction  of  solutions  to  problems  including  non¬ 
linear  interactions  of  electromagnetic  waves  with  an  ionized  atmosphere.  The  generation  of 


3-9 


electro-acoustic  modes  may  also  be  considered  as  a  phenomenon  contributing  to  power 
absorption  of  a  transmitting  antenna  in  a  reentry  environment. 

3.3  ANTENNA  BREAKDOWN 

This  section  discusses  theoretical  methods  of  predicting  the  electric  field  intensity  at  which 
microwave  breakdown  of  the  atmosphere  takes  place.  Breakdown  electric  field  is  a  useful 
parameter  because  it  can  be  related  to  a  minimum  number  of  independent  variables.  On  the 
other  hand,  the  power  handling  capability  of  an  antenna  depends  on  the  field  distribution  of 
the  particular  antenna. 

The  discussion  covers  static  breakdown  and  reentry  effects.  Static  breakdown  is  defined  as 
dielectric  breakdown  In  an  undisturbed  environment.  Ttr  subject  of  reentry  effects  is 
concerned  with  breakdown  of  an  antenna  mounted  on  a  vehicle  in  hypersonic  flight  in  the 
atmosphere. 

3.3.1  STATIC  BREAKDOWN 

The  subject  of  microwave  breakdown  in  gases  has  been  thoroughly  surveyed  by  MacDonald  in 
52 

a  recent  book.  This  book  presents  practically  all  of  the  basic  theory  of  high  frequency  gas 
discharges.  Moreover,  it  indicates  how  breakdown  field  strengths  can  be  calculated  reliably. 
The  book  is  thus  useful  in  providing  the  basic  physical  theories  on  the  subject  as  well  as  the 
practical  tools  needed  for  applications.  In  addition,  a  whole  chapter  is  devoted  to 
experimental  methods. 

Let  us  now  give  a  brief  summary  of  basic  breakdown  theory,  beginning  with  a  qualitative 
description  of  the  breakdown  process.  The  impressed  electromagnetic  field  exerts  force  on 
any  free  electrons  which  may  be  present  in  the  gas.  Through  the  process  of  collisions  of 
electrons  with  neutral  gas  particles,  the  field  does  work  on  the  electrons,  raising  their 
average  temperature  well  above  that  of  the  gas.  Ionization  of  neutral  particules  by  impact 
of  high  energy  electrons  competes  with  electron  attachment  to  neutral  particles  to  give  a  net 
rate  of  electron  population  gain.  This  net  ionization  rate  competes  with  losses  by  diffusion 
until  the  point  is  reached  where  the  ionization  builds  up  at  a  very  rapid  rate  and  the  gas 
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becomes  conducting.  The  various  collision  processes  are  very  important  in  determining 
breakdown,  in  that  they  control  the  energy  distribution  as  well  as  the  number  density  of 
electrons. 


The  kinetic  theory  approach  to  breakdown  prediction  treats  ,  he  electron  distribution  function 
F  (v,  r,  t),  which  represents  the  number  of  electrons  with  velocities  close  to  v,  in  a  small 
volume  of  space  at  the  point  r,  and  at  the  time  t.  This  distribution  f  action  obeys  the 
Boltzmann  equation53’ 54 


57  +  v  .  ?F  +  a  •  V  F  =  C 
ot  v 


(3.  3.1-1) 


where  V  is  the  gradient  in  configuration  space,  ^  is  the  gradient  in  velocity  space,  a  is  the 
acceleration,  and  C  represents  the  effects  of  collisions. 

The  collision  term  C  may  be  written  as  an  integral  involving  the  distribution  function  and  the 

54 

collision  cross  section,  and  hence  knowledge  of  collision  cross  section  as  a  function  of 
electron  velocity  is  sufficient  in  principle  to  solve  the  Boltzmann  equation.  However,  the 
difficulties  involved  are  so  great  that  accurate  results  are  possible  in  only  a  few  cases. 

Hydrogen  and  helium  have  practically  constant  collision  frequencies  above  an  energy  of  about 
3  or  4  electron  volts,  and  this  fact  has  been  used  by  MacDonald  and  Brown55’  56  and  Reder 
and  Brown  to  calculate  breakdown  fields  for  these  two  gases  and  for  admixtures  of  mercury. 
The  results  agree  very  well  with  experiment,  especially  for  the  case  of  helium  containing  a 
small  amount  of  mercury.  The  presence  of  mercury  serves  to  convert  the  metastable 
excitation  level  of  helium  to  ioniz-  t  on  of  mercury  by  an  exchange  reaction.  This  has  the 
effect  that  the  gas  has  no  excitation  levels  below  the  effective  ionization  level,  greatly 
simplifying  the  analysis.  The  theoretical  results  for  neon  and  neon-argon  mixtures,  for 
which  the  collision  frequency  is  approximately  proportional  to  the  square  root  of  the  electron 
energy,  do  not  agree  with  experiment  quite  so  well  as  for  constant  collision  frequency  gases. 


Most  practical  microwave  breakdown  theory  is  based  on  phenomenological  models,  because 
of  the  difficulties  encountered  with  collisional  processes  in  the  kinetic  theory  approach.  In 
the  phenomenological  theory,  the  electron  continuity  equation  is  written  in  terms  of  macro¬ 
scopic  rate  constants,  which  can  be  determined  experimentally  as  functions  of  electric  field 
ard  pressure.  The  electron  continuity  equation  is  written  as 


-♦ 
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(3.  3. 1-2) 
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where  n  is  the  electron  density,  T  is  the  electron  flux,  and  P  is  the  net  effect  of  all  sources 
and  sinks. 


Free  diffusion  ;>?  electrons  applies  to  static  breakdown,  where  the  initial  electron  density  is 
very  low.  Then 


r  -  -  V(Dn)  (3.  3. 1-3) 

where  D  is  the  diffusion  coefficient  When  the  net  ionization  rate  per  electron  is  represented 
by  the  symbol  v,  the  continuity  equation  may  be  written  as 
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The  Townsend  criterion,  carried  over  from  the  theory  of  direct  current  discharges,  has  been 
shown  by  Herlin  and  BrownJ  to  apply  to  microwave  discharges.  According  to  this  criterion, 
the  breakdown  point  is  reached  when  the  right  hand  side  of  Equation  3.  3. 1-6  is  zero,  since  a 
slightly  higher  electric  field  would  cause  the  runaway  production  of  electrons. 

1 

i 

When  either  the  electric  field  or  the  electron  mobility  is  not  uniform  in  space,  as  would  be 

the  case  in  antenna  breakdown,  the  diffusion  equation  must  be  solved  in  the  form  of  Equation 

3.  3. 1-4,  with  the  left  hand  side  set  equal  to  zero.  However,  it  is  possible  to  use  an  effective 

diffusion  length  which  can  be  derived  from  the  solution  of  the  diffusion  equation.  For  example, 

H  *rlin  and  Brown  derived  the  solution  for  a  TM  -mode  cylindrical  cavity  and  MacDonald 
6  0  ^ 

and  Brown  derived  the  results  for  a  spherical  cavity. 

52 

It  can  be  shown,  through  certain  considerations  in  the  kinetic  theory,  that  the  important 
pai ameter  with  respect  to  electric  field  E  and  the  microwave  radian  frequency  ^  is 


(3.  3.1-7) 

where  Ee  is  the  effective  field  and  v q  is  the  collision  frequency.  The  effective  field  concept 

is  useful  in  the  phenomenological  theory  of  breakdown,  as  in  the  analysis  of  Brown  61  Rose 

62  63  ' 
and  Brown,  and  Gould  and  Roberts.  The  Townsend  criterion  is  written  as 

L_ 

"i  ~  "a  +  a2  (3.  3. 1-8) 

A 

I 

where  and  v ^  are  the  ionization  and  attachment  frequencies,  respectively.  The  first 
Townsend  coefficient  a  is  used  in  place  of  according  to  the  equation 


v.  =  a  u  E 
i  e 


(3. 3. 1-9) 
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where  /t  is  the  dc  mobility.  The  ratio  of  diffusion  coefficient  to  mobility  is  invoked  to  give  the 
relation 


o'  _  B  ,  2  u 

p  p  3  (pA)  (E  A)  (3#  3- 1_ 

e 

where  p  is  pressure,  u  is  average  electron  energy,  and  /3  is  the  number  of  attachments  per 

centimeter  of  path  length.  Data  exists  for  a/p,  /3/p,  and  u  as  functions  of  E/p  in  air.^’^ 

Hence,  Equation  3.  3. 1-10  may  be  used  to  give  Ee/p  as  a  function  of  pA,  provided  E  can  be 

e 

calculated. 


The  difficulty  connected  with  the  phenomenological  theory  is  that  the  effective  field  is  not  a 

61 

sound  concept,  since  the  collision  frequency  is  not  constant.  Brown  obviates  this  difficulty 
by  calculating  an  effective  constant  collision  frequency  from  the  measured  ac  mobility,  using 


(e/m)  v 

c 


(3.  3. 1-li) 


where  e/m  is  the  ratio  of  charge  to  mass  for  the  electron.  He  gets  a  value  of  u  equal  to 
9  c 

4.  3  x  10  p.  The  results  of  breakdown  calculations  then  agree  with  experimental  data  fairly 

well  over  a  limited  range  of  pressures. 

It  appears  that  the  most  reliable  and  generally  useful  method  of  calculating  breakdown  is  that 
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given  by  MacDonald.  This  method  is  a  phenomenological  one  which  uses  the  effective  field 
concept  in  a  way  which  is  much  less  sensitive  to  the  specific  form  of  collision  frequency 
assumed  than  is  the  above  method.  The  diffusion  coefficient  can  be  calculated  from  the 
variation  of  collision  frequency  with  electron  energy,  if  the  electron  energy  distribution  is 
assumed;  however,  it  can  be  shown  that  the  form  of  the  distribution  function  has  little  effect 

/*  r 

on  the  result.  Hence  the  approximation  uf  MacDonald,  Gaskell,  and  Gitterman  works  well. 
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(3.  3. 1-12) 


/  0.  9  E  \ 

Dp  =  (29  +  — -  e  1  io4  cm2  -  Torr/sec 

This  equation  was  derived  from  assuming  a  Unear  dependence  of  Dp  on  average  electron 
energy  in  terms  of  the  latter's  experimentaUy  determined  dependence  on  E/p.  Equation 
3.  3. 1-12  is  equivalent  to 


where  X  is  the  wavelength  and 


S 


pX 


29  + 


0.  9  EX 

(pX)  +  (35.  6)2 


where  EX  is  in  volts  and  pX  is  in  Torr-cm. 


(3.  3. 1-14) 


The  net  ionization  rate  v  may  be  related  to  the  electric  field  and  the  pressure  by  using 
experimental  breakdown  data.  The  result  is  shown  in  Figure  8-18  of  Reference  52  in 
terms  of  v  X  versus  pX  for  fixed  values  of  EX.  Figure  8-19  in  Reference  52  shows  a  plot 
of  Equation  3. 3. 1-14.  The  solution  to  a  given  problem  in  continuous  wave  breakdown  in  air 
is  found  by  the  points  of  intersection  of  these  two  sets  of  curves,  as  shown  in  Figure  8-20  in 
Reference  52,  since  at  breakdown,  pX  is  equal  to  DX/A  2.  The  curves  are  superimposed 
such  that  the  point  S  =  1  coincides  with  the  point  j/X  =  104  (X/A  )2.  Thus,  when  the  wavelength 
and  the  diffusion  length  are  known,  the  breakdown  field  versus  pressure  can  be  determined 
quite  easily  by  this  method. 

The  above  method  of  calculating  continuous  wave  breakdown  can  be  generalized  to  include 
pulse  (or  finite  time)  effects,  following  the  approach  of  Gould  and  Roberts. 63  For  example, 

for  pulsed  breakdown  at  low  pulse  repetition  rate,  the  net  ionization  frequency  is  determined 
by 
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UbT  ~  in  <Vno)  (3.3.1-15) 

where  Tis  the  pulse  length,  is  the  electron  density  required  for  breakdown,  and  n  is  the 
initial  electron  density.  The  value  of  n^  is  usually  taken  as  the  critical  value,  1013/X2,  where 
X  is  in  cm  and  n^  is  in  cm  ;  and  n  is  a  fairly  small  number.  For  a  given  pulse  length, 
Equation  3.  3. 1-15  determines  p  ,  which  is  actually  related  to  v  of  the  CW  analysis  by 

D 

*b  ~  V  ~  .2  (3.  3. 1-16) 

A 

The  procedure  is  to  estimate  E  from  continuous  wave  theory  and  to  determine  D  from  Equation 
3.  3. 1-13.  Then  v  is  determined  from  Equation  3.  3. 1-16.  This  gives  E  from  CW  theory,  and 
the  process  can  be  iterated  to  determine  E  more  accurately.  In  case  of  rapid  pulse 
repetition  rates,  the  decay  time  of  the  afterglow  between  pulses  must  be  taken  into  account. 
Then  the  breakdown  criterion  is  that  the  buildup  of  ionization  during  the  pulse  is  just  equal  to 
the  decay  between  pulses.  A  slightly  greater  net  ionization  rate  during  the  pulse  would  give  a 
runaway  condition.  As  the  pulse  repetition  rate  is  increased,  the  breakdown  field  approaches 
the  continuous  wave  case. 

3.3.2  REENTRY  EFFECTS 

The  effects  of  reentry  on  the  breakdown  conditions  are  to  change  the  rates  at  which  the 
various  breakdown  controlling  processes  proceed  and  to  add  processes  not  considered  in  the 
static  breakdown  case.  The  following  paragraphs  will  discuss  some  of  the  changes  mentioned 
in  the  literature  and  will  give  an  estimate  of  their  importance  in  affecting  breakdown  where 
present  information  allows. 

An  important  effect  of  the  reentry  induced  hypersonic  flow  field  is  an  apparent  shift  in 
altitude  due  to  increased  air  density52.  It  is  important  to  note  that  the  air  density  is  the 
proper  independent  variable  in  determining  breakdown,  rather  than  pressure,  since  the  air 
is  relatively  warm. 
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Another  effect  of  reentry  is  the  change  in  gas  composition.  Due  to  the  high  levels  of  gas 
temperature  and  pressure  resulting  from  the  imparting  oi  vehicle  kinetic  energy  to  its 
environment,  the  thermodynamic  gas  state  lor  the  reentry  environment  is  significantly 
different  from  the  room  temperature  state.  This  results  in  a  shift  in  the  population  densities 
from  and  (principal  constituents  of  room  temperature  air)  to  the  dissociated  forms 
1,0  and  O  (which  occur  in  appreciable  concentrations  in  a  hypersonic  cone  boundary  layer). 
The  resultant  lowering  of  impact  ionization  potentials, which  may  be  seen  from  the  following 
brief  table, will  increase  the  ionization  rate,  enhancing  breakdown. 


Species 

ev 

°K  x  10“ 

NO 

9.25 

10.7 

°2 

12.2 

14.2 

N2 

15.6 

18.1 

N 

14.5 

10.8 

O 

13.6 

15.8 

An  additional  effect  of  the  high  temperature  environment  may  be  a  substantial  increase  in 

66 

attachment  rate  for  the  reentry  environment.  Loss  of  electrons  by  attachment  to 

.  67 

molecular  oxygen  has  been  found  to  be  most  important  relative  to  other  constituents  of  air. 

66 

Epstein  notes  that  the  cross  section  for  this  process  may  be  significantly  changed  when 
oxygen  is  vibrationally  excited  at  high  temperature,  resulting  in  a  substantial  increase  in  the 
attachment  rate  for  the  reentxy  environment. 

The  second  main  difference  in  the  reentry  environment  is  the  presence  of  a  significant  level 
of  ambient  ionization.  This  increased  ion  population  derives  from  three  effects: 

a.  Intense  heating  in  the  bow  shock  wave  and  viscous  boundary  layer  about  the  vehicle 
during  hypersonic  flight 

b.  Introduction  of  easily  ionizable  impurities  in  the  flow  field  due  to  heat  shield  ablation 

c.  Ambient  ionization  of  gases  in  the  upper  atmosphere. 
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These  phenomena  force  reconsideration  of  two  aspects  of  the  breakdown  analysis.  The  first 
aspect  of  the  problem  is  simply  the  possibility  of  increasing  the  initial  plasma  frequency  by 
virtue  of  the  higher  electron  density  such  that  it  becomes  an  appreciable  fraction  of  the 
applied  signal  frequency.  This  condition  tends  to  shorten  the  time  required  to  reach  critical 
electron  density  and  therefore  breakdown.  This  factor  must  then  be  taken  into  account  by 
adding  a  term  to  the  electron  continuity  equation  representing  ionization  rate  in  the  absence  of 
an  electromagnetic  field.  This  term  is  determined  by  the  flow  field  calculations  which 
consider  the  thermal,  ambient,  and  contaminant  ionization  rates  as  a  function  of  gas  state. 

The  second  aspect  of  the  existence  of  appreciable  positive  ion  densities  is  to  cause  the 

66 

diffusion  process  to  become  ambipolar  rather  than  free  as  in  the  classical  case.  Allis  and 
68 

Rose  have  analyzed  the  transition  from  free  to  ambipolar  diffusion.  They  indicate  that  the 

transition  occurs  when  the  Debye  length  for  electrons  becomes  smaller  than  the  diffusion 

52 

length.  Since  ambipolar  diffusion  is  slower  than  free  diffusion  by  a  factor  of  about  forty  , 
the  rate  of  loss  of  free  electrons  will  be  reduced,  resulting  in  a  significant  reduction  in  power 
required  for  breakdown. 

Other  terms  which  may  be  required  for  a  description  of  reentry  induced  breakdown  are  the 

66  69 

rates  of  recombination  and  detachment.  According  to  Lin  and  Tear  the  dissociative 

recombination  of  electrons  and  nitric  oxide  ions  (the  most  abundant  ion  present  in  high 

—3  “*3/2  3  66 

temperature  air)  proceeds  at  a  rate  of  approximately  3  x  10  T  cm  /sec.  Epstein 

G 

shows  that  although  according  to  the  data  of  Muschlitz,  detachment  is  negligible  in  the  region 
from  3000°  to  5000°  K,  further  work  is  required  before  a  definite  conclusion  can  be  made 
regarding  the  importance  of  detachment. 

Another  effect  of  reentry  is  the  additional  electron  loss  mechanism  caused  by  the  convection 

70 

of  air  past  the  antenna.  Kelly  and  Margenau  introduce  this  difference  by  imting  that  even 

for  continuous  wave  breakdown  the  atmosphere  is  only  subject  to  microwave  heating  for  a  time 

T  =  L/V  (L  is  the  length  of  the  *>^cenna  and  V  is  the  velocity  of  the  antenna).  Thus  breakdown 

can  occur  only  if  the  electron  concentration  becomes  large  enough  in  time  T  since  the  antenna 

71 

encounters  a  new  environment  every  T  seconds.  Fante,  in  order  to  take  account  of  this 
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effect,  couples  the  velocity  of  the  gas  flowing  past  the  antenna  V  with  the  electron  continuity 

70 

equation  by  rewriting  the  term  dn/dt  as  V  •  Vn.  Kelly  and  Margenau  also  redefine  break¬ 
down  for  a  moving  medium  as  that  point  in  time  when  the  signal' frequency  w  equals  the 

66 

plasma  frequency  w  .  Epstein  points  out  that  this  definition  is  probably  superior  (though 
P 

it  is  only  approximately  true)  to  the  classical  one  where  the  power  required  for  breakdown 
is  defined  as  the  incident  power  level  which  causes  the  net  ionization  rate  to  become  positive. 
This  classical  definition  fails  to  predict  changes  in  the  electron  density  if  convection  is 
present,  and  thus  provides  no  information  regarding  the  way  in  which  the  electron  density 
approaches  its  final  value. 


Finally,  the  effect  of  nonuniform  flow  field  distributions  has  been  considered  by  Fante  *  in 
terms  of  its  effect  on  the  diffusion  equation.  This  effect  is  taken  into  account  by  solving  the 
diffusion  equation  for  the  given  distribution  of  the  diffusion  coefficient. 


It  is  probably  best  to  calculate  the  effect  s  of  reentry  on  antenna  breakdown  by  the  pheno¬ 
menological  approach,  rather  than  the  kinetic  theory  approach.  Classical  breakdown  in  air, 
as  discussed  above,  is  successfully  predicted  on  this  basis.  It  should  then  be  possible  to 
formulate  the  approach  in  such  a  way  as  to  obey  a  correspondence  principle  with  respect  to 
the  classical  limit.  In  other  words,  the  results  should  approach  the  classical  results  when 
the  reentry  conditions  become  insignificant.  This  is  probably  most  conveniently  done  by 
rewriting  the  continuity  equation  with  additional  terms  representing  each  of  the  reentry 
effects,  or  with  modified  terms  in  certain  cases. 


3.4  NONHNEAR  ATTENUATION  AND  DISTORTION 

This  section  deals  with  the  problem  of  interaction  between  a  partially  ionized  plasma  medium 
and  incident  high  power  radio-frequency  electromagnetic  radiation.  The  interaction  induces 
a  change  in  the  equation  of  state  of  the  plasma  due  to  field  energy  coupling  into  the  electron 
gas.  Asa  result,  the  electron  distribution  function  is  altered  and  becomes  dependent  on  the 
electric  field  intensity.  Consequently,  transport  coefficients,  such  as  thermal  and  electrical 
conductivities,  which  are  deduced  from  the  distribution  function,  become  dependent  on  the 
local  field  intensities.  While  various  forms  of  nonlinearities  may  be  exhibited  from  the 


interaction  of  electromagnetic  radiation  with  a  plasma,  the  particular  class  covered  in  this 

section  is  restricted  to  the  case  where  the  incident  field  deposits  energy  in  the  plasma  and 

modifies  the  electron  temperature,  electron  density,  and  effective  collision  frequency.  The 

present  interest  is  in  the  propagation  of  electromagnetic  radiation  in  a  plasma,  and  the 

conductivity  coefficient  is  of  paramount  importance.  More  specifically,  if  f  and  f  represent 

e  i 

the  electron  and  ion  velocity  distribution  functions,  then  the  components  of  the  current  i  are 

k 

related  to  the  electric  field  components  E^  by 

In  the  above  equation  v^  is  the  k  component  of  the  velocity  vector  and  o  Is  the  scalar 
conductivity  coefficient.  In  the  presence  of  a  magnetic  field  the  scalar  conductivity  is 
replaced  by  a  tensor  conductivity  coefficient.  For  most  practical  purposes  the  ion  conduction 
current  may  be  ignored  and  then  Equation  3. 4-1  is  approximated  by, 


j 


k 


,3 

v,  d  v 
k 


a  E,. 
k 


(3.  4-2) 


In  the  preceding  equation  the  subscript  on  the  electron  distribution  function  has  been  dropped 
for  convenience.  If  the  distribution  function  is  a  nonlinear  function  of  the  electric  field 
intensity,  the  conductivity 

0  _  i  -effvki\ 

Et  ’  Ek  <3'4"3 


becomes  a  function  of  the  electric  field  intensity. 

In  the  following  literature  survey  the  electromagnetic-plasma  nonlinear  interactions  v/ill  be 
grouped  into  three  parts  depending  on  the  geometrical  configuration  of  the  plasma  medium, 
namely,  infinite,  semi-infinite,  and  slab. 
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3.  4. 1  INFINITE  PLASMA  MEDIUM 

The  propagation  of  electromagnetic  radiation  in  an  unbounded  nonlinear  plasma  medium 

excludes  boundary  effects  and  allows  emphasis  to  be  placed  on  the  physics  of  the  interaction 

process.  An  extensive  investigation  of  this  phase  of  the  phenomena  has  been  carried  out  by 

72  J 

Ginzburg  and  Gurevich.  In  fact,  their  paper  provides  an  appropriate  introduction  and 

background  to  the  entire  field.  They  show  that  the  nonlinear  interaction  of  electromagnetic 

waves  with  partially  ionized  plasmas  can  be  characterized  in  terms  of  several  basic  relaxation 

times.  These  are  as  follows: 

T  =  Electron  density  relaxation  time 

Ten  =  Electron  temperature  relaxation  time 

^em=  27r/(j  =  Period  of  incident  electromagnetic  radiation 

^ee  ~  Relaxation  time  for  the  distribution  function  to  become  Maxwellian 
2  , 

=  L  /D  =  Characteristic  diffusion  time 


L  =  Characteristic  diffusion  length 

D  =  Diffusion  coefficient 

From  the  relative  order  of  magnitude  of  the  relaxation  times  it  is  possible  to  infer  the  specific 
nonlinear  response  of  a  plasma  to  intense  electromagnetic  radiation.  Thus  if 

(3.4-4) 
(3. 4-5) 


T  »  T  , 
em 


T  »  T  , 
en  em 
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the  electron  density  and  temperature  cannot  follow  the  rapid  variations  of  the  electromagnetic 
field  but  assume  a  steady  state  value  proportional  to  the  mean  square  time  average  of  the 
electric  field.  On  the  other  hand  if 


T  «  t 

em 


(3. 4-6) 


r  «  t 
en  em 


(3.  4-7) 


then  variations  in  electron  density  and  temperature  have  ac  components  similar  to  those  of  the 

electric  field  intensity.  This  phenomenon  exhibits  the  well  known  "Luxembourg"  effect.  A 

further  subdivision  of  the  response  takes  place  depending  on  the  relaxation  time  T  ,  or  the 

ee 

time  required  for  the  electrons  to  reach  a  Maxwellian  distrioution.  Thus  two  possibilities 
exist  depending  on  whether  T  is  greater  or  less  than  the  time  required  for  the  relaxation  of 

vv 

e.octron  temperature.  A  further  subdivision  depends  on  whether  the  electron  density  ci  anges 
or  is  constant  under  the  influence  of  the  external  field. 
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Ginzburg  and  Gurevich  follow  Allis  and  obtain  a  solution  to  the  Boltzmann  equation  by 
expanding  the  distribution  function  in  a  spherical  harmonic  series  or 

V  vi 

f  ■  fo+  w  + - '  (3.4-8) 

where 

f  =  Isotropic  part  of  the  distribution  function 
.  th 

f.  =  i  component  of  the  first  spherical  component  of  the  distribution  function 
th 

Vj  =  i  component  of  electron  velocity 
w  =  Electron  speed 


3-22 


It  is  found  that  when  the  thermal  velocity  is  greater  than  the  drift  velocity,  then  it  is 
sufficient  to  maintain  only  two  terms  in  the  series.  It  is  advantageous  to  use  the  expansion 
m  spherical  harmonics  rather  than  the  perturbation  method,  since  with  the  former  technique 
it  is  possible  to  provide  an  accurate  treatment  of  the  energy  dependence  of  the  electron- 
neutral  and  electron-ion  collisions  for  momentum  transfer.  By  taking  the  steady-state  case 
@  *o^t  =  0)  and  the  following  collisions, 

a.  Elastic  electron-neutral 

b.  Elastic  electron-ion 

c.  Inelastic-neutral 

d.  E  lectr  or  -electron 

Ginzburg  and  Gurevich  arrive  at  two  expressions  for  the  isotropic  distribution  function  f  for 
the  electrons:  °  • 


a-  If  Tee  <<C  Ten  then  iQ  assumed  a  Maxwellian  form. 


b*  U  ree  >;>  Ten’  111611  electr°n-electron  collisions  do  not  affect  the  form  of  f  and 


w 


fQ  =  C  exp  /- 


^  f  kT  + 


mw  dw 


e  E  /3m  G^2  +  y  (w)l 


(3. 1-9) 


Furthermore,  when 


e2  E2/3m  G  Co,2  +  y2  <w)]  «  kT,  4_ 

Equation  3. 4-9  reduces  to  a  Maxwellian  form.  In  the  preceding  equations  G  is  the  relative 
fractional  loss  of  energy  of  an  electron  per  inelastic  collision  (assumed  small  such  as  in 


excitation  of  rotational  or  vibrational  states)  and  T  is  the  gas  temperature.  Equation  3. 4-9 
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was  originally  derived  by  Margenau  . 

The  ac  conductivity  is  obtained  from  the  first  spherical  harmonic  component  of  the  electron 
distribution  function  f^',  where 


w  (w)  -  "en  <w)  +  "ei  <">.  (3.4-14) 

where  and  represent  velocity  dependent  electron-neutral  and  electron-ion  collisions. 

Thus  from  Equation  3. 4-13  if  fQ  is  a  function  of  the  electric  field  intensity  then  o,  and 
therefore  the  propagation  constant,  is  a  function  of  E.  The  appropriate  value  for  f  is  given 
by  Equation  3. 4-9  or  the  Maxwellian  form  at  temperature  T  , 
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One  of  the  earlier  approaches  to  the  nonlinear  interaction  problem  was  by  Margenai  ,  He 
considered  an  infinite  medium  and  arrived  at  the  distribution  function  (isotropic)  gi  ven  by 
Equation  3. 4-9  or 


f 

o 


w 


Cexp 


•/ 


mw  dw 


2  2  2  2 
kT  +  e  E  /3m  G  Tco  +u  (w)D 


(3.  4-15) 


Here  T  was  taken  greater  than  T  and  interelectron  collisions  were  not  considered.  The 
ee  en 

electron  density  was  constant, and  elastic-ion  and  inelastic  electron-neutral  collisions  were 
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ignored.  Epstein  has  also  examined  a  similar  problem. 
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Sodha  and  Palumbo  studied  the  propagation  of  an  ellipitically  polarized  wave  in  a  nonlinear 
homogenous,  and  unbounded  magnetoplasma.  The  electron  distribution  function  was  not 
Maxwellian.  Electron-ion  and  electron-el  jctron  collisions  were  not  considered.  Electron 
density  was  constant.  The  components  of  the  conductivity  tensor  were  found  to  be  functions 
of  the  square  of  the  local  electric  field  amplitude. 
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Sodha  and  Palumbo  have  also  examined  the  problem  where  several  electromagnetic  waves 
propagate  in  an  unbounded  partially  ionized  plasma.  The  isotropic  part  of  the  electron 
distribution  function  was  derived  and  found  to  be  dependent  on  the  frequency  and  square  of  the 
electric  field  amplitude  for  each  wave.  In  the  limit,  for  a  single  wave,  the  isotropic 
distribution  function  reduced  to  the  case  where  only  one  wave  is  present.  The  conductivity 
depends  on  the  frequency  and  square  of  the  electric  field  amplitude  for  each  wave  present  in 
the  plasma.  Thus  the  propagation  characteristics  of  any  single  wave  depend  on  the  power 
density  of  the  other  waves  propagating  in  the  medium. 
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The  problem  of  multi-frequency  propagation  has  been  considered  by  H,  Mott  in  conjunction 
with  antenna  breakdown  in  air.  The  breakdown  depends  not  only  on  frequency  but  also  on  the 
envelope  c'  the  electric  field.  The  envelope  may  be  determined  by  phasor  addition  of  the 
individual  components.  A  minimum  level  is  obtained  for  n  equal  amplitude  signals  by  in-phase 
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addition.  The  resulting  envelope  peaks  determine  breakdown  by  acting  as  a  continuous  wave 
breakdown  field,  or  if  E  ff  is  the  effective  dc  field, 


Eeff  "  n  Eo  "  'W  (3- 4-16) 

where  Eq  is  the  electric  field  intensity  of  each  signal.  The  total  minimum  power  for  break¬ 
down  is  given  by, 

P/Peff  =  n  (ft)  =  1/n  (3.4-17) 

'  eff ' 

Additional  work  on  multi-frequency  signals  has  been  carried  out  by  Fante  and  Mullin79.  Here 
an  attempt  is  made  to  determine  the  effective  electric  field  strength  for  a  composite  signal 
such  as  noise.  For  a  stationary  random  signal  E  (t),  the  effective  field  is  given  by 


(3. 4-18) 


where  is  the  velocity  independent  collision  frequency,  and  S  (f)  is  the  spectral  density. 

The  transmission  of  power  by  the  noise  signal  is  compared  to  that  of  a  single  frequency 
component.  At  high  collision  frequencies  there  is  no  difference  between  the  signals.  However, 
when  -  o>0  the  single  frequency  can  transmit  more  power  than  the  noise  signal.  Here 
w0  =  +  w2)/2,  where  ^  and  w  represent  the  upper  and  lower  cutoff  frequencies. 

3. 4. 2  SEMI-INFINITE  PLASMA 

The  interaction  of  electromagnetic  radiation  with  a  semi-. infinite  plasma  introduces  boundary 
phenomena  and  facilitates  the  analytical  investigation  of  the  behavior  of  spatial  inhomogeneities 
in  the  plasma.  This  problem  was  considered  by  King80,  Yen81,  and  Papa82’ 83.  The  result 
of  the  perturbation  (electric  field)  modifies  the  boundary  between  the  nonionized  and  ionized 


medium  in  such  a  way  that  it  appears  to  move  towards  the  radio  frequency  source.  The 
electron  temperature  and  electron  density  become  functions  of  both  time  and  space.  King 
and  Yen  make  use  of  the  time-independent  WKB  approximation.  They  assume  that  the 
isotropic  part  of  the  electron  velocity  distribution  function  is  Maxwellian,  but  do  not 
account  for  the  velocity  dependence  of  the  electron-neutral  cross  sections.  No  heat  conduction 
in  the  electron  gas  is  assumed.  Papa's  analysis  is  based  on  the  three  moments  of  the 
Boltzmann  aquation  for  the  electron  gas.  The  moment  equations  are  truncated  by  setting  the 

heat  flow  vector  proportional  to  the  temperature  gradient.  The  relaxation  time  is  short 
enough  for  f  to  remain  Maxwellian.  Ionization  by  electron  impact  on  neutrals,  electron- 

electron  attachment,  and  electron-ion  recombination  are  included  in  the  analysis.  The  con¬ 
ductivity  is  given  by  Equation  3.4-13.  The  continuity  and  energy  balance  equations  are 
solved  numerically  for  the  electron  density  N  and  temperature  T  .  It  is  also  shown  by  Papa 

G  © 

that  the  time  scale  T  is  much  shorter  than  T. 
en 
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Dysart  applied  the  theoretical  approach  given  by  Papa  *  and  obtained  numerical  solutions 

to  the  time  dependent  problem.  His  results  almost  completely  duplicate  Papa's  original 
presentation.  Extensive  graphical  data  is  presented. 
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Epstein  treated  the  semi-infinite  problem  for  a  chemically  reacting  and  inhomogeneous 
plasma  medium.  He  did  not  take  into  account  electron-ion  collisions,  and  the  velocity 
dependence  of  electron-neutral  collisions,  and  vhe  form  of  the  distribution  function.  He 
assumed  that  the  ratio  of  electromagnetic  energy  transferred  to  the  plasma  to  the  internal 
energy  (shock  induced)  was  small  and  thereby  linearized  the  fluid  dynamical  and  chemical 
equations.  Under  these  somewhat  restricted  conditions  he  showed  that  a  small  electro¬ 
magnetic  field  produced  substantial  changes  in  the  reflection  coefficients. 
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Papa  has  also  considered  the  case  where  a  dc  magnetic  field  is  imposed  parallel  to  the 
electron  density  gradient.  Time  dependence  was  taken  into  account.  He  showed  that 
gradients  in  electron  density  cause  the  incident  field  to  experience  frequency  shift.  The 
frequency  shift  of  the  left-hand  circularly  polarized  wave  differs  from  the  frequency  shift  of 
the  right-hand  wave.  Dysart8^  has  also  supplemented  this  work  by  performing  numerical 
computations. 


3.4.3  PLASMA  SLAB 

The  nonlinear  interaction  between  electromagnetic  radiation  and  an  inhomogeneous  plasma 

slab  represents  a  problem  which  in  many  ways  closely  simulates  practical  conditions.  Papa 
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and  Case  have  examined  the  steady-state  nonlinear  response  of  an  inhomogeneous  plasma 

slab  when  illuminated  by  intense  plane  wave  radiation.  The  constituents  were  selected  to 

correspond  to  those  of  air  behind  a  strong  shock.  The  constituents  included  are  the  neutrals 

N  »  NO,  N,  and  O.  It  is  shown  that  for  a  large  range  of  shock  parameters  T  «  r  , 

*  ee  en 

T  <K  T  ,and  f  is  Maxwellian,  Thus  the  rate  coefficients  can  be  defined  in  terms  of  an 
ee  D  o 

electron  temperature  T  .  The  energy  balance  equation  is  solved  numerically  and  the  results 

e  2  2 

plotted  graphically  for  T  versus  E  /u  .  It  is  shown  that  it  is  possible  to  have  multivalued 

®  2  2 

solutions  of  T  for  certain  ranges  of  E  /w  .  In  this  range  the  electron  gas  can  exist  in  three 
6 

temperature  states.  The  reflection  and  transmission  coefficients  were  computed  by  Papa 

and  Case  employing  the  Runge-Kutta  technique  to  integrate  Maxwell’s  equations.  The 

2  2 

reflection  coefficient  is  presented  as  a  function  of  E,  for  several  values  of 

inc  “-pfmax) 

Here  E^  is  the  incident  field  and  wp^maxj  is  the  maximum  plasma  frequency  of  the  initial 
electron  density  profile.  It  is  found  that  the  reflection  coefficient  is  relatively  constant  until 
a  critical  value  is  reached  corresponding  to  the  electron  avalanche.  At  this  point  the 
reflection  coefficient  increases  sharply  while  the  transmission  coefficient  drops  to  low  values. 
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Dysart  has  obtained  numerical  solutions  for  electromagnetic  transmission  through  a  plasma 
slab  under  steady  state  conditions.  The  plasma  is  taken  as  inhomogeneous  and  isotropic  with 
parameters  corresponding  to  the  flow  around  a  blunt  vehicle  travelling  at  18, 000  ft/sec  at  an 
altitude  of  200, 000  feet.  The  constttutents  are  Nr  NO,  and  O  at  5000°K.  The  radio  frequency 
is  taken  much  greater  than  collision  frequencies  so  that  the  conductivity  given  by  Equation 
3. 4-13  reduces  to, 


a 


2  2 
N  e  v  N  e 
+1_e_ 

*  moj 
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(3.4-19) 


and 
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(3.  4-20) 
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with 


v.  (T  )  +  u, 
i  e  ion 


(T  ), 

e 


(3. 4-21) 


where  u.  (Tg)  is  the  effective  collision  frequency  for  the  component  of  the  neutral  constituent. 
Nume  deal  data  is  presented  graphically  and  includes  reflection  and  transmission  coefficients. 

The  effort  by  Dysart  is  strictly  computational  and  the  theory  is  identical  to  that  presented 
by  Papa  and  Case8'. 


The  case  of  propagation  through  an  inhomogeneous,  nonlinear  magnetoactive  plasma  slab  for 
time  invariant  electron  density  was  examined  by  Papa88.  The  components  of  the  conductivity 
tensor  are  functions  of  the  power  density  of  three  modes,  namely,  right-hand,  left-hand,  and 
longitudinal.  For  normal  incidence  an  elliptically  polarized  incident  field  launches  only  a 
right-hand  and  s  left-hand  mode  into  the  plasma  for  a  dc  magnetic  field  parallel  to  the  electron 
density  gradient  These  two  modes  are  uncoupled  at  low  power.  At  higher  input  power,  the 
modes  do  not  propagate  independently  of  one  another.  If  desired,  it  is  possible  to  compute  the 
net  reflection  and  transmission  coefficients  for  the  left-  and  right-hand  modes. 

3.5  ALLEVIATION  TECHNIQUES 

Some  methods  of  alleviating  reentry  communicationc  blackout  are  briefly  surveyed,  including 
aerodynamic  shaping,  reducing  the  free  electron  density,  and  altering  the  propagation 
properties  of  the  plasma. 

The  general  principle  in  aerodynamic  shaping  is  either  to  locate  the  antenna  outside  of  the 
plasma  sheath  or  to  modify  the  flow  field  so  that  the  plasma  density  is  reduced  in  the  vicinity 


of  the  antenna.  Three  techniques  are  compared  in  Reference  89.  The  communications  fin  is 
a  narrow  angle  wedge  fin  projecting  beyond  the  bow  shock  of  the  vehicle.  Since  the  fin  is 
narrow,  the  plasma  generated  by  it  is  weak  compared  to  the  plasma  generated  by  a  blunt 
vehicle.  Thus  if  the  antenna  system  were  located  on  the  fin,  the  effects  of  plasma  attenuation 
could  be  alleviated.  The  problems  involved  include  structural  stability  of  the  fin,  and 
aerodynamic  stability  of  the  vehicle  with  such  a  fin. 

The  aerospike  antenna  is  a  narrow  antenna  projecting  into  the  plasma  sheath.  A  coolant  is 
injected  into  the  flow  stream  upstream  of  the  antenna. 

Flow  diverters  are  small  projections  from  the  vehicle  into  the  flow  stream  which  would 
modify  the  flow  stream  so  that  the  sheath  thickness  over  the  area  of  the  antenna  is  reduced. 
Flow  diverters  would  be  lighter  than  the  two  previously  mentioned  techniques,  and  would 
require  less  modification  to  existing  vehicle  designs. 

Reference  89  also  considers  plasma  alleviation  by  coolant  injection.  Among  various  materials 
considered,  water  was  chosen  as  most  effective.  The  jet  of  water  would  be  shot  normal  to 
the  surface  of  the  vehicle.  A  shock  would  form  because  of  the  interaction  of  the  jet  an i  the 
flow  field,  and  the  jet  would  break  up  into  droplets.  The  water  would  then  cool  the  plasma 
by  evaporation. 
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Ohio  State's  Electro-Scicnce  Laboratory  is  designing  a  system  which  would  project  a  jet  of 
cool  cas  ahead  of  the  vehicle.  Not  only  would  this  system  cool  the  plasma  formed,  but  it  would 
also  protect  the  vehicle  since  the  shock  would  form  on  the  tip  of  this  jet,  rather  than  on  the 
nose  of  the  vehicle.  The  gas  spike  is  to  be  tested  on  a  trailblazer  vehicle. 

The  injection  of  an  electronegative  gas  such  as  SF  could  have  three  desirable  effects  as 

o 

follows:  (1)  it  would  act  as  a  coolant  because  of  the  heat  of  dissociation  of  SF  into  its 

o 

component  atoms.  (2)  It  could  reduce  the  equilibrium  concentration  of  free  electrons  due  to 
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attachment.  (Good  points  out  that  the  predominant  mechanism  is  actually  a  three-body 
reaction  where  the  third  body  carries  off  the  excess  association  energy. )  (3)  It  serves  as  a 


catalyst  for  air  recombination.  Thus  it  can  be  used  in  conjunction  with  another  coolant  to 
hasten  equilibrium. 

Reference  89  also  considered  injecting  dust  into  the  plasma,  so  that  the  dust  particles  could 
serve  as  nuclei  for  electrons  to  collect  on.  One  dust  particle  could  clean  out  a  volume  with 
a  radius  of  one  Debye  length.  Since  the  smallest  particle  which  could  be  easily  handled  was 
thought  to  be  on  the  order  of  1  micron,  it  was  concluded  that  other  methods  would  weigh  less. 

Attenuation  for  right-hand  circularly  polarized  waves  propagating  normal  to  a  plasma  layer 

will  be  reduced  in  the  presence  of  a  magnetic  field  in  the  direction  of  propagation.  This  is 

because  the  Lorentz  force  checks  the  ability  of  the  electron  generated  field  to  cancel  the 
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imposed  field.  Rothman  and  Morita  consider  the  problem  analytically  for  rectangular 
coordinates.  Hie  conditions  for  enhancement  are  that  the  electron  cyclotron  frequency  be 
greater  than  the  imposed  angular  frequency,  and  that  the  square  of  the  collision  frequency  be 
much  less  than  the  square  of  the  difference  between  the  cyclotron  frequency  and  the  imposed 
angular  frequency.  Rothman  and  Morita  carried  out  experimental  laboratory  work  using  an 
RF  discharge  which  demonstrated  qualitative  enhancement  and  polarization  in  the  presence  of 
a  magnetic  field. 

93  94 

Samaddar  *  considers  magnetic  window  effects  in  special  cases  of  nonrectangular 

geometries.  In  particular,  he  considers  cylindrical  and  conical  geometries  which  approximate 

vehicle  shapes.  His  method  utilized  calculating  the  dielectric  tensor  as  a  function  of  the 

applied  magnetic  field,  and  the  plasma  and  collision  frequencies.  To  gain  effective  control 

over  attenuation,  mutual  orientation  of  the  antenna  and  imposed  magnetic  field  must  be  chosen 

so  that  the  field  components  are  independent  of  components  of  the  dielectric  tensor  parallel  to 
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the  imposed  magnetic  field.  This  approach  has  been  analyzed  in  terms  of  weight  require¬ 
ments  and  found  to  be  feasible  only  if  a  greater  system  loss  figure  than  with  the  other  alleviants 
(20  db  rather  than  10  db)  could  be  tolerated.  The  use  of  superconducting  magnets  hag  some 
promise  for  short  flights,  but  cooling  equipment  becomes  a  limiting  factor. 
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A  critical  analysis  of  these  and  other  techniques  should  be  made  before  any  firm  conclusions 
can  ne  drawn.  Fluid  mechanic  effects  of  shaping  and  additives  will  be  considered  more 
systematically  in  future  reports.  Alleviation  of  antenna  breakdown  will  be  given  more 
emphasis  nlso. 
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SECTION  4 
CONCLUSION 


It  has  been  determined  that  the  construction  of  an  accurate  nonequilibrium  flow  field  is 
possible  through  the  judicious  use  of  the  nonequilibrium  streamtube  solution  in  conjunction 
with  a  perfect  gas  or  equilibrium  flow  field  solution.  The  method  is  approximate  only  to  the 
extent  that  the  bow  shock  and  streamwise  pressure  distributions  are  approximate. 

The  problem  of  analyzing  antenna  near-field  distributions  cannot  be  set  up  in  a  perfectly 
general  way  which  will  give  good  results  in  a  variety  of  cases.  However,  the  general 
approaches  used  in  the  theory  of  the  linear  antenna  provide  a  good  basis  for  approximate 
approaches  to  more  general  problems. 

A  good  deal  of  applicable  literature  exists  on  the  linear  effects  of  plasma  on  antenna  impedance 
and  pattern  distortion.  Somt'  of  these  theories  could  prove  useful  in  making  parametric 
calculations. 

The  prediction  of  antenna  breakdown  by  the  kinetic  theory  abroach  is  not  feasible  except  for 
certain  ideal  cases.  Phenomenological  theories  have  succeeded  in  predicting  the  results  of 
experiment  for  static  breakdown  in  air.  Whether  or  not  the  phenomenological  theory  can  be 
extended  to  cover  the  effects  of  reentry  remains  to  be  seen. 

Nonlinear  plasma  effects  and  some  techniques  for  alleviation  of  linear  plasma  effects  have 
been  surveyed.  Assessment  of  the  importance  of  these  two  technical  areas  with  respect  to 
the  present  application  must  await  further  evaluation. 
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i  ’3.  » E-ST  R  ACT 

This  report  is  the  first  quarterly  technical  report  under  this  study  of  reentry  effects 

6  ^an8mlS8lon  of  electronic  countermeasure  (ECM)  signals.  Illustrative 
results  for  tile  aerothermochemical  properties  are  presented  for  typical  slender-body 

1  ™TtryA  “n8‘  Th6Se  re8UltS  WCre  cbtained  from  311  approximate  nonequilibrium  ~  I 
inviscid  shock  layer  flow  field  solution,  which  utilizes  a  combination  of  reacting  streamlube 
and  cTuiL.brium  flow  field  solutions.  A  discussion  of  the  calculation  method  is  included.  I 

in  thG  area  °f  tbBOry  °f  antGDnaS  311(1  Propagation  are  summarized.  ’  . 

S  field  d^ribbZr  T  <Xlmm0n  USeare  1118011886(1  in  Volu^eIL)  Methods  of  calculating 
near  field  distributions  of  antennas  are  discussed  in  a  general  way.  Theoretical  analyses  of 

linear  plasma  effect,  are  surveyed.  Antenna  breakdown  is  considered  in  terms  of  botii  the 

stax  c  air  environment  and  the  reentry  environment.  The  theories  of  nonlinear  propagation 

effects  are  summarized  and  finally  a  review  of  some  alleviation  techniques  is  given. 

In  general,  it  appears  that  nearly  all  of  the  important  reentry  effects  on  antennas  and 
propoagation  can  be  reasonably  well  characterized  using  existing  theory.  However 
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